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In this letter, I have proposed a model of variable modified Chaplygin gas and shown its role in 
accelerating phase of the universe. I have shown that the equation of state of this model is valid 
from the radiation era to quiessence model. The graphical representations of statefinder parameters 
characterize different phase of evolution of the universe. All results presented in the letter concerns 
the case = 0. 



PACS numbers: 

Recent observations of the luminosity of type la supernovae indicate [1, 2] an accelerated expansion of 
the universe and lead to the search for a new type of matter which violate the strong energy condition 
p + 3p < 0. The matter consent responsible for such a condition to be satisfied at a certain stage of 
evolution of the universe is referred to as dark energy. There are different candidates to play the role 
of the dark energy. The type of dark energy represented by a scalar field is often called Quintessence. 
In particular one can try another type of dark energy - the so-called pure Chaplygin gas which obeys an 
equation of state like [3] 

P^-B/p, (1) 

where p and p are respectively the pressure and energy density and i? is a positive constant. Subse- 
quently the above equation was modified to the form (known as generalized Chaplygin gas) 

p = -B/p" with < a < 1. (2) 

This generalized model has been studied previously [4, 5]. There are some works on modified Chaplygin 
gas obeying an equation of state [6, 7] 

p — Ap with < a < 1, A, B are positive constants. (3) 

This equation of state shows radiation era (when A = 1/3) at one extreme (when the scale factor a{t) 
is vanishingly small) while a ACDM model at the other extreme (when the scale factor a{t) is infinitely 
large). Guo and Jhang [8] first proposed variable Chaplygin gas model with equation of state (1), where 
i3 is a positive function of the cosmological scale factor 'a' i.e., B = B{a). This assumption is reasonable 
since B{a) is related to the scalar potential if we take the Chaplygin gas as a Born-Infeld scalar field [9]. 
Later there are some works on variable Chaplygin gas model [10]. 

The metric of a homogeneous and isotropic universe in FRW model is 

dr^ 



ds^ ^dt^ ~a\t) 



+ r'^{de'^ + sin^edcf^) 



1 — fcr-^ 

where a{t) is the scale factor and k (= 0,±1) is the constant curvature of their spatial sections. 
The Einstein field equations are 



(4) 



fc 1 



(5) 
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and 

^ = > + 3p) (6) 

where p and p are energy density and isotropic pressure respectively (choosing SttG = c = 1). 
The energy conservation equation is 

/9 + 3^(/9 + p) = (7) 

Now, I have introduced variable modified Chaplygin gas with equation of state (3) where is a positive 
function of the cosmological scale factor 'a' (i.e., B = B{a)) as 

p = Ap- with < a < 1, ^ is constant > 0. (8) 

At all stages it shows a mixture. Also in between there is also one stage when the pressure vanishes 
and the matter content is equivalent to a pure dust. 

Now, for simplicity, assume B{a) is of the form 

B{a) = Boa-"" (9) 
where Bq > and n are constants. Using equations (7), (8) and (9), I have the solution of p as 



3(l + a)So 1 C 



{3{1 + a){l + A) - n} a" a3(i+^)(i+") 



(10) 



where C > is an arbitrary integration constant and 3(1 + A)(l + a) > n, for positivity of first term. 
Here n must be positive, because otherwise, a ^ oo implies p ^ oo, which is not the case for expanding 
Universe. 

Prom equation (5), for A; = 0, we get the explicit form of t in terms of a as 

^1 C 
t = KaW+^) 2Fi[ -z,l-z,~— a-^(TO] 

z(i + a) K 

where 

n 



(l+a)" 



/ n 




GBqz 





2{1 + a){3{l + A){1 + a) - n} 
The deceleration parameter q has the expression 

a 



For accelerating universe, q must be negative i.e., a > i.e., 

2{l + a)-n ^3(i+„)(i+^)_„ C{1 + 3A) 
3(l + a)(l + A)-n Wq 

which requires n < 2(1 + a). Since < a < 1, this implies n < 4. 



(11) 



This expression shows that for small value of scale factor we have decelerating universe while for large 
values of scale factor we have accelerating universe and the transition occurs when the scale factor has 
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C(l+3A){3{l+a){l+A)-n} 
3Bo{2(l+a)-n} 



the expression a = 

Now for small value of scale factor a{t), I have 



3(l+a)(l + A)-n 



^ c— 

P- a3(i+A) ' 

which is very large and corresponds to the universe dominated by an equation of state p = Ap. 
Also for large value of scale factor a{t), 



(12) 



i{l + a)Bo 



{3(l + a)(l + A)-n} 



a i+o 



and I obtain 



p= -1 + 



(13) 



(14) 



3(1 + a), 

which correspond to quiessence model (i.e., dark energy with constant equation of state) [13]. 

Note that n = corresponds to the original modified Chaplygin gas scenario [7] , in which the modified 
Chaplygin gas behaves initially radiation and later as a cosmological constant. However, equation 

(10) shows that, in the variable modified Chaplygin gas scenario, it interpolates between a radiation 
dominated phase {A = 1/3) and a quiessence-dominated phase described by the constant equation of 



state p = jp where 7 = — 1 



3(l+a) 



I have described this particular cosmological model from the field theoretical point of view by intro- 
ducing a scalar field (p and a self-interacting potential V{(l)) with the effective Lagrangian 



1 :r 



-V{cP) 



(15) 



In the paper of Gorini et al [4], it has been shown that for the simple fiat Friedmann model with 
Chaplygin gas can equivalently described in terms of a homogeneous minimally coupled scalar field (j). 
Following Barrow [11], Kamenshchik et al [3, 12] have obtained homogeneous scalar field <j){t) and a 
potential F(^) to describe Chaplygin cosmology. 

Now, I consider the energy density and pressure corresponding to a scalar field <j) having a self- 
interacting potential V^(^). The analogous energy density p^ and pressure for the scalar field are 
the following: 



</>2 + 1/(0) =p = 



3(1 + a)Bo 



{3{1 + a){l + A) - n} 



(16) 



and 



V{c^)=Ap--^ 

pa 



3(l + a)Bo 



1 

— + 



C 



+ a){l + A) - n} a" a3(i+A)(i+a) 



- -Boa 



3(1 + a)Bo 



C 



mi + a){l + A) - n} a" a3(i+>i)(i+a)_ 



(17) 



Hence for flat universe (i.e., = 0), I have 
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^/T+A 



{3(l + a)(l + ^)-n} 



, , ./i / (^/x{u + y) + ^/y{u + x)) 

2 log(V^+ - ^ log V ___ 1 



and 



(18) 



V{cP) = -{l-A) 



3(1 + a)Bo 



C 



{3(l + a)(l+A)-n} a" a3(i+^)(i+") 



3(l + a)So 



+ 



C 



{3(l + a)(l + A) -n} a" o3(i+^)(i+«) 



(19) 



where x = j^, y = 3{l + a) and u = ^ (| - l) 



7"(l-i). 



The graphical representations of (f> against a and V{(f>) against (f> have been shown in figures 1, 2 and 
figures 3, 4 respectively for A= 1/3. Figures 1 and 3 show the fixed value of a = 0.6 and various values 
of n (= 1, 2, 3). In this case, V{<p) increases with (f) and slope of the curves decreases as n increases. 
Figures 2 and 4 show the fixed value of n = 2 and various values of a (= 0.9, O.5., 0.3). In this case also, 
y(</)) increases with </> and slope of the curves decreases as a decreases. Figures 1 and 2 describe the 
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Fig. 5 shows the variation of s against r for different values of n (= 3, 1) and for a = 0.6, A = 1/3. 



scalar field always decreases with the evolution of the universe. 

In the paper [4], the flat Friedmann model filled with Chaplygin fluid has been analyzed in terms 
of the recently proposed ^'statefinder'^ parameters [14]. The statefinder diagnostic along with future 
SNAP observations may perhaps be used to discriminate between different dark energy models. The 
above statefinder diagnostic pair i.e., {r,s} parameters are constructed from the scale factor a{t) and its 
derivatives upto the third order as follows: 



where H is the Hubble parameter and q (= — fr) is the deceleration parameter. These parameters 
are dimensionless and allow us to characterize the properties of dark energy in a model independent 
manner. The statefinder is dimensionless and is constructed from the scale factor of the Universe and its 
time derivatives only. The parameter r forms the next step in the hierarchy of geometrical cosmological 
parameters after H and q. 

For Friedmann model with flat universe (i.e., k = 0), 

4 = (21) 



and 



a 1 3p 



So from equation (20) we get 



Thus, I get the ratio between p and p: 
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I = fcil (24) 

For variable modified Chaplygin gas, using equations (7), (23) and (24), I get the relation between r 
and s: 



18(r - l)s^ + lSas{r - 1) + 4a(r - 1)^ = 9sA(l + a){2r + 9s - 2) + 3ns(2r - 2 - QsA) (25) 

Figure 5 shows the variation of s with the variation of r iov A = 1/3 and for a = 0.6 and n = 3, 1. 
The portion of the curve on the positive side of s which is physically admissible is only the values of 
r greater than {l + +^)}- The part of the curve between r = 1 and r = 1 + |A(1 + A) with 
positive value of s is not admissible (we have not shown that part in the figure 5) because for the Chaply- 
gin gas under consideration we face a situation, where the magnitude of the constant Bq becomes negative. 

Thus the curve in the positive side of s starts from radiation era and goes asymptotically to the dust 
model. But the portion in the negative side of s represents the evolution from dust state (s = — oo) to 
the quicsscncc model. Thus the total curve represents the evolution of the universe starting from the 

radiation era to the quiessence model. 

In this work, I have presented a model for variable modified Chaplygin gas. In this model, I am 
able to describe the universe from the radiation era [A = 1/3 and p is very large) to quiessence model 
(p is small constant). So compare to Chaplygin gas model, the present model describe universe to a 
large extent. Also if, I put A = O with a = 1, then I can recover the results of the Chaplygin gas 
model. If put n = then variable modified Chaplygin gas model reduces to only modified Chaplygin 
gas model [7]. In figure 5, for {r,s\ diagram the portion of the curve for s > between r = 1 to 
r = 1 + 1^(1 + A) is not describable by the modified Chaplygin gas under consideration. For example, if 
I choose r = 1.03, A = 1/3 then from the curve s = 0.01 which corresponds to g = 3/2 and hence I have 
from the equation of state, _B < which is not valid for the specific Chaplygin gas model considered 
here. At the large value of the scale factor I must have some stage where the pressure becomes negative 
and hence B has to be chosen positive. It follows therefore that a portion of the curve as mentioned 
above should not remain valid. 
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